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Abstract 

Determination of the nuclear parton distributions within the framework of perturbative 
QCD, the DGLAP equations in particular, is discussed. Scale and flavour dependent 
nuclear effects in the parton distributions are compared with the scale and flavour 
independent parametrizations of HIJING and of the Hard Probe Collaboration. A 
comparison with the data from deep inelastic lepton-nucleus scattering and the Drell- 
Yan process in proton-nucleus collisions is shown. 
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1 Introduction 



In a high-energy collision of hadrons or nuclei A and B, inclusive cross sections for the 
production of a particle c involving a large scale Q ^> Aqcd, can in the leading twist 
approximation be computed by using collinear factorization, 



da(Q 2 , y/s) A B^c+x = 



Z A f! /A (x u Q 2 ) + (A-Z A )fr /A (x u Q^ 

hj=q,q,g 

Z B ff B (x 2 , Q 2 ) + (B- Z B )ff B (x 2 , Q 2 )\ ® da(Q 2 , x u x 2 )^ c+x (1) 



where da(Q 2 , x±, x 2 )ij^ c+x is the perturbatively calculable differential cross section for 
the production of c at the scale Q, X\ )2 ~ Q/y/s are the fractional momenta of the 
colliding partons % and j, and /f is the distribution of parton flavour i in a proton 
(neutron) of the nucleus A, and correspondingly fp n '' B is that for the nucleus B. The 
number of protons in A(B) is denoted by Z A (Z B ). In the leading twist approximation, 
on which we shall focus in the following, multiple scattering of the bound nucleons does 
occur but all collisions are independent of each other and only one-parton densities 
are needed. At this level all the possible higher twist terms, suppressed by 1/Q 2 
but enhanced by the nuclear geometry (thickness of the nuclei), are neglected. As 
correlations between partons are not considered, the nuclear effects enter only via 
the nuclear parton distribution functions (nPDF) fP> n / A > B , The nPDF differ from 
the parton distributions of the free proton but obey the same DGLAP [|l] evolution 
equations. The DGLAP evolution of the nPDF has been studied in e.g. 0-0. For 
studies of next-to-leading twist factorization involving two-parton distributions, see 



10, 11 



The cross sections of hard processes measured in deeply inelastic I A collisions and in 
pA collisions offer the experimental constraints necessary for pinning down the parton 
densities in nuclei. Once the nPDF are known, the reference cross sections for hard 
probes of dense matter in ultrarelativistic nucleus-nucleus collisions can be computed 
and used in the search of the signals of the QGP. Another motivation for the study 
of the nPDF is that the measurements of certain hard processes in nuclear collisions 
are useful also for extracting information of the parton distributions of the free proton, 



provided that the contribution of the nuclear effects can first be removed [12]. 

The purpose of this note is to discuss the constraints of the nuclear parton densities 
within the DGLAP framework, such as the ones in [7|, ]8|, We shall compare the 
outcome of the studies 0, ||, the EKS98 parametrization of the nuclear effects, with 



scale and flavour independent parametrizations of HIJING |I3| and the Hard Probe 
Collaboration [14], and especially, with the data. 
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1.1 DIS and nPDF 



The cleanest way of getting information of the nPDF is from deeply inelastic lepton- 
nucleus scattering (DIS) experiments. The differential cross section for deeply inelastic 
Ip scattering in the one-photon exchange approximation can be expressed as 

W ff (*, q») f _ + 1 1 + tftftf 1 (2) 



dQ 2 dx Q 4 x I * a ' 2 l + i^Q 2 ) 

with the standard Lorentz-invariant variables x = Q 2 / (2p • g) and y = p ■ q/p ■ k, where 
p, g and k are the four-momentum of the proton, the exchanged virtual photon, and 
the incoming lepton, correspondingly Virtuality of the photon is Q 2 = — q 2 , mass of 
the nucleon is M, and E\ is the energy of the incoming lepton in the target rest frame. 
In principle < x < A but for the discussion here the small cumulative tails of the 
distributions at x > 1 can be safely neglected. 

From the measurements with different nuclear targets, one knows that the ratio 
of the absorption cross sections of longitudinal and virtual photons with the target 
nucleons, R(x, Q 2 ) = al /a? does not significantly depend on the target nucleus A at 



scales Q > 1.5 GeV |TR ITBf- The ratio of cross sections with different targets thus 



reflects the ratio of the nuclear structure functions FA/F B : 

da lA /dQ 2 dx F 2 A (x, Q 2 ) Z A F v 2 IA {x, Q 2 ) + (A - Z A )F? /A (x, Q 2 ) 



da^/dQ 2 dx F 2 B (x, Q 2 ) Z B F? /B (x, Q 2 ) + (B - Z B )F? /B (x, Q 2 



(3) 



where p^ ,n ' A,B are the structure functions of bound nucleons. The ratio of the DIS 
cross sections from IA and ID is then related to the average structure functions per 
nucleon as 

E A (T Q2) _ \da^/dQ 2 dx ^ \FA IjFT + F? /A ) + |(f - 1)(F? A - F^ A ) 
Fal M ] ~ \da lD /dQ 2 dx ~ \F? iffl/D + F n/ D) > W 

where the numerator is written as a sum of isospin symmetric and non-symmetric 
terms. The nuclear effects in deuterium are small, less than 1 percent, so to a first 
approximation these can be neglected. Measurements of the ratio Rp 2 (x,Q 2 ) have 
revealed clear and systematic nuclear effects in different regions of Bjorken-x ||17||-|28|: 



"shadowing"; a depletion at x<0.1, 
"anti-shadowing" ; an excess at 0.1 < x < 0.3, 
"EMC effect"; a depletion at 0.3 < x < 0.7, 
"Fermi motion" ; an excess towards x — > 1 and beyond. 
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Systematics of Rp 2 in A and in x has been extensively studied e.g. in the experiments 
EMC [0, 0, U , SLAC S H, BCDMS [||, NMC E665 |7|, ||. Since the 

Q 2 dependence of Rp 2 is quite weak at x > 0.1 it has been more difficult to probe. Data 
with high enough precision, however, exist: the NMC has some years ago discovered 



a clear Q 2 -dependence in the ratio of the cross sections da^ n /da^ c [P6fl , i.e. the scale 
dependence of the ratio F 2 Sn /F 2 c , at x ~ 0.01. 



2 The DGLAP analysis 

Generally, perturbative QCD (pQCD) cannot predict the absolute parton distributions. 
However, once the starting distributions are given at a scale Qo, pQCD successfully 
predicts the evolution in Q 2 in the form of the DGLAP evolution equations 0]. The 
global DGLAP analyses of the parton distributions of the free proton, such as MRS 
[29] and CTEQ fl30"| |, involve a determination of those input distributions which, when 



evolved to different (higher) values of Q, give the best overall agreement with the data 
from different hard processes. Conservation of momentum and baryon number are 
maintained by the DGLAP evolution, and they are used as further constraints. 

In the leading twist framework the situation is exactly the same for the nPDF: the 
data from hard processes in nuclear collisions at various values of x and Q, together 
with the momentum and baryon number conservation, can be used to constrain the 
nonperturbative input distributions of partons in bound protons at some initial scale 
Qo 3> Aqcd- The link between the scales Qo and Q is given by the DGLAP equations 
in the whole range of x studied. 

In general, the perturbative QCD scale evolution of the nPDF has been extensively 
studied in the literature, see e.g. 0-@ and |JT], |32|. The origin of the nuclear effects 



is an interesting question but is beyond the scope of the DGLAP analysis. It has been 
suggested that for the DGLAP evolution one may compute the initial conditions at a 
scale Qo from a model, as is done e.g. in || ||, and then apply the DGLAP equations 
to describe the evolution in Q. Nevertheless, even in this case the key feature is the 
detailed comparison with the existing data. DGLAP analyses of the nPDF which 
attempt to rely only on the data in the determination of the initial conditions, are 
presented in Refs. ]7|, Q. 

2.1 Quarks and antiquarks 

In the QCD- improved parton model (in leading order, or in the DIS-scheme in any 
higher order) the structure function F 2 of the proton (neutron) can be written in terms 
of its parton distributions as 



if " )/A (z, Q 2 ) = £ e 2 xf^ A (x, Q 2 ) + xft )/A (x, Q 2 ) 

q=u,d,s,... 



(5) 
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As in the case of the free nucleons, the parton distributions of bound neutrons in 
isoscalar nuclei are obtained through isospin symmetry, /"^ = f^t and f^t = 
This is expected to be a good approximation for non-isoscalar nuclei as well. 

It is convenient to define the nPDF for each parton flavor i through the modifications 
of the corresponding distributions in the free proton, 

R?^Q 2 ) = f O X n*\ (6) 
n (x, Q 2 ) 

where we shall assume that the parton distributions of the free proton are fully known. 
For example, below the mass threshold of the charm quark, we can write 

da / ^ _ 5Qa + u A + d A + d A ) + As A + (^f - 1)3Qa + u A - d A - d A ) 

where u A = f '^ A = R A (x,Q 2 )f^(x,Q 2 ), and similarly for the other quarks. The ratio 
Rp 2 measured at these scales thus constrains the individual ratios R A and R A in certain 
combination. 

Obviously, more constraints are needed in order to pin down the ratios Rf. These 
can be obtained from the measurements of the Drell-Yan dileptons in pA collisions 
by E772 j33] and E866 |34j] which offer high enough precision for getting statistically 
significant constraints. In the lowest order, the ratio of the differential cross sections 
for the Drell-Yan process in pA and pD collisions is given by 

A , 2 s _ \da p DY /dx 2 dQ 2 



\do\^ Y I 'dx 2 dQ 2 

{A[u x {u A + d A ) + u^v* + d A )\ + [d^di + u A ) + d x {d A + u A )\ + As lS A + ...}/N DY 
+(— - l){A[ Ul (u A - d A ) + u x {u A - d A )\ + [d x {d A - u A ) + d x (d A - 4)}}/N DY (8) 

where the invariant mass of the lepton pair is Q 2 , and the subscript 2 (1) refers to the 
fractional momentum x 2 (xi) of the parton from the target (proton). The denominator 
is 

N DY = 4[u x (u 2 + d 2 ) + u x (u 2 + d 2 )\ + [d x (d 2 + u 2 ) + d x (d 2 + u 2 )\ + As x s 2 + ... (9) 

with the dots denoting the heavier flavours. Again, the ratios R A and R A are probed 
in certain combination. The typical x and Q 2 range probed by the measurements of 
the ratios Rp 2 and Rjj Y can be seen in Fig. 1. 
In addition, conservation of baryon number, 



3 = f A dx £ \fv/ A (x,Q 2 )-f? /A (x,Q* 



q=u,d 
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Figure 1: The correlation between the Bjorken x and the photon virtuality Q 2 in deep 
inelastic [iA scattering measured by the NMC |23|, 24, [25|] and E665 |27], ||] (denoted by 
DIS). The same (with x = x<i) for the measurements of the Drell-Yan process in pA by E772 



(denoted by DY). The starting scale of the DGLAP analysis 



is Ql 



Jl dx £ Q 2 )f^ Q 2 ) - Rf(^ Q 2 )f%(x, Q 2 ) 



q=u,d 



(10) 



can be used to pin down the valence quark distributions || ||. On the r.h.s. of Eq. 
(|10|), the small cumulative tails at x > 1 have been neglected. 

Even in an ideal case, where the experimentally measured DIS and DY ratios of Eqs. 
(|) and (|J) would lie along a constant scale Qo in a wide range of x, the number of 
experimental constraints above would not be enough to fully fix the ratios Rf(x,Ql) 
and Rq(x,Q 2 ). In reality, the situations is even more difficult: the data on both DIS 
and DY are given only in distinct regions, within which the values of x are strongly 
correlated with Q 2 , as illustrated in Fig. |l|. Therefore, a recursive procedure similar 
to that in the global analyses of the parton distributions of the free proton has to be 
adopted in order to determine the initial ratios Rf(x,Ql). Notice that in comparison 
with the free proton case, an additional variable, the mass number A, appears. 

To simplify the determination of the input nuclear effects for valence and sea quarks 
(without invoking any specific model) one may in a leading approximation assume 
them to be separately flavor- independent: R^ v (x,Ql) « R^ v (x,Ql) ~ Ry(x,Ql), and 
R£(x,Ql) w R$(x,Ql) w Rf(x,Ql) w R s (x, Ql) g], 0. Note that this approximation 
is needed only at Ql but the observation in is that it remains good also in the 
evolution to higher Q 2 . In this approximation the problem reduces to constraining the 
three independent input ratios, Ry, Rj; and Rq at the initial scale Qq. The details 
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of an analysis using this approach can be found in M, below we only summarize the 
available constraints in each region of x. 



At x> 0.3 the valence quark distributions dominate the structure function F 2 A , 
and Rp 2 ~ R v . The DIS data for Rp therefore only constrains the magnitude 
of the EMC effect and the Fermi-motion in R v but not in R$ or in R£. For the 
sea quarks, it is assumed in |7J that Rs(x > 0.3, Qq) ~ Ry(x > 0.3, Qq) (cf. the 
discussion for gluons below). To what extent the Drell-Yan production measur- 
able in nuclear collisions at the SPS could probe the EMC effect of the nuclear 
sea quarks, was recently studied in Ref. f35fl. 



At 0.04 < x < 0.3 the DIS and DY data both constrain the ratios R$ and R v but 
from different regions of Q 2 , as shown in Fig. [[]. In addition, Ry is restricted 
by conservation of baryon number. An outcome of the DGLAP analyses P, |], 0] 
is that at the input scale Qq ~ 2 GeV 2 the sea quark content of a nucleus 
remains smaller than that of the free proton, i.e. that no antishadowing appears 
in Rj(x,Q 2 ). 

At x < 0.04 the DIS data for the ratio i?4 extend down to x ~ 5 ■ 10 -3 in the 
region Q > 1 GeV relevant for the DGLAP analysis. In the analyses ||, [7j, |D| the 
nuclear valence quarks have less shadowing than the sea quarks, which is mainly 
due to the conservation of baryon number. Also more strongly shadowed valence 
quarks have been suggested ||. 

The DIS data for Rp 2 at x < 5 • 10~ 3 only exist in the region Q < 1 GeV which 
can be considered nonperturbative and not treatable with the DGLAP equations 
(see Fig. [l]). A saturation behaviour, flattening of Rf 2 in x — > 0, is observed 
along the experimentally probed values of Q 2 |24], ^ . Such a behaviour, a weak 
dependence of Rp 2 on x, can also be expected at Qq, provided that the sign of 
the slope of the Q 2 dependence of Rp 2 in the nonperturbative region remains the 



same (positive) as what is measured at x ~ 0.01 in the perturbative region p6 
The DIS data in the non-perturbative region should thus give a lower bound for 
Rp 2 (x, Qq) at small values of x. The sea quarks dominate over the valence quarks 
in this region, so practically only the ratio Rg is constrained by the DIS data. 

2.2 Gluons 

For the gluons the situation is less straightforward since direct measurements of the 
gluon distributions in nuclei are difficult. In principle the measurements of D-mesons 
in pA collisions at various cms-energies can be used for pinning down the nuclear gluon 
distributions pp] , [37]. These measurements already exist but so far the error bars in the 
results of experiment E789 [|38| are too large for getting a stringent constraint for Rq. 
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The future measurements of D production in pA collisions by the NA60 experiment 
at SPS, by PHENIX at RHIC and hopefully also by ALICE or CMS at the LHC, 
will provide very important input for fixing the gluon distributions in nuclei J^J . Also 
direct photons in nuclear collisions can be used for this purpose, possibly also diffractive 
scattering in DIS [39]. Production of J/^fr in pA always involves strong final state effects 



(J/ty suppression in normal nuclear matter), which makes the extraction of the initial 
state effects from the data very difficult 

In the DGLAP analysis, the best (but still indirect) constraint on the ratio Rq(x, Q 2 ) 
is provided by the measurements by NMC |26| of the Q 2 dependence of the ratio 
F 2 Sn /F 2 c . The scale evolution of F 2 is coupled to the gluon distributions at small values 



of x, where gluons dominate, approximately as [fiT 



dF 2 (x,Q 2 )/dlogQ 2 » ^ xg (2x,Q 2 ). (11) 
Based on Eq. (|TT|), the Q 2 slope of Rp 2 at small values of x becomes 

dR U x ,Q 2 ) „ 10a s xg(2x,Q 2 ) { A 2 A 2 A 
d\ogQ 2 ~ 277T F 2 d (x,Q 2 )\ Kg[2X ' Q } H ^ Q >}' W 

which suggests that the more deeply gluons are shadowed, the slower is the evolution 
of Rp 2 . So far only the NMC data p6f of the Q 2 dependence of the ratio F 2 Sn /F 2 c has 



sufficient precision for getting stringent constraints, as first pointed out in Ref. |42| . 
Yet another indirect constraint, the momentum sum rule, 

l = J A dx £ xff /A (x,Q 2 )^ fdx Y, Rf(x,Q 2 )xf[(x,Q 2 ) (13) 

^ i=g,u,u,... i=g,u,u,... 

should be taken into account in the DGLAP approach. A few percent flow of mo- 
mentum from quarks and antiquarks to gluons is expected relative to the free pro- 
ton case 0, IJ. Since Jq 1 cLxRq(x, Ql)xg p (x, Ql) > dx xg p (x,Qf ) ), antishadowing 
Rq(x, Qq) > 1 is bound to exist in some region of x. Below we summarize the con- 
straints on Rq(x, Ql) in different regions of x. 

• At 0.02 < x < 0.2 the NMC data on the Q 2 evolution of the ratio F 2 Sn /F 2 c |(J set 
the main constraint through the DGLAP evolution (the full equations; Eq. ([T2|) 
holds at the small values of x only). These NMC data extend down to x — 0.0125, 
so only gluons at x>0.02 can be constrained. Especially, it is observed that the 
Q 2 slope of the ratio Fl ? n / F 2 C is clearly positive at small values of x, indicating 
that obviously also the Q 2 slope of Rp 2 (x,Q 2 ) is positive. Then, according to 
Eq. (|12 ), Rq(2x,Q 2 ) > Rp 2 (x,Q 2 ). The data thus seems to rule out the case 
where the shadowing of nuclear gluons is much stronger than the shadowing 
observed in Rp- In a solution consistent with the data on F 2 Sn /F 2 c and with the 
momentum sum rule the gluons have less shadowing than the sea quarks and 
more antishadowing than the valence quarks or in Rp 2 0, [|2 . 
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At x< 0.02 stringent experimental constraints do not exist for the gluons at the 
moment. Assuming, however, that the Q 2 slope of Rp 2 remains positive, the 
measured saturation of shadowing in Rp at Q 2 <C 1 GeV 2 [^|, [27| gives the 
lower bound for the shadowing of Rp 2 (x, Qq). The weak x dependence of Rp 
(together with dRp 2 /dQ 2 > 0) indicates also a weak x dependence of Rq(x, Qq), 
and it is then concievable to expect that Rq(x,Qq) ~ Rp 2 (x,Ql) for x <C 1 0. 
The evolution is stable in the sense that the approximate equality remains true 
within about 5 % even after the DGLAP evolution from Q ~ 1 GeV to Q ~ 100 
GeV ||. In order to pin down the gluons in this region, high-precision DIS 
measurements at small values of x (but Q 2 > 1 GeV 2 ) would be needed. More 
constraints in this region of x are expected from the measurements of dileptons 
originating from the decays of the D mesons in PHENIX at RHIC [RBI E7| and 



in ALICE or CMS at the LHC [[37]]. 



At x > 0.2 there are currently no clear experimental constraints available for the 
gluons. Conservation of momentum does not reveal whether an EMC effect exists 
for the gluons or not: about 30 % of the gluon momentum comes from x>0.2, 
so, say, a 10 % net change in the momentum content of the EMC region can 
be compensated by roughly a 6% net effect in the region where the antishadow- 
ing bump is anticipated and which contains about half of the gluon momentum. 
Thus, the amount of antishadowing is not affected to the extent that it would vi- 
olate the constraints obtained from x < 0.2 (see the estimates of the uncertainties 
in [0). In the DGLAP evolution equations, the valence quarks act as source of 
gluons in the gluon evolution, and the gluons in turn feed the sea quark evolu- 
tion. In course of the scale evolution, the experimentally verified EMC effect of 
valence quarks will generate a similar EMC effect for the gluons which in turn 
transmit the effect into the sea quark distributions, as seen in H, where no input 
EMC effect was assumed for the gluons. A consistent assumption in the DGLAP 
framework therefore is to include an EMC effect already for the initial ratios 
Rq and Rg. In this way the nuclear modifications Rf remain stable against the 
evolution, i.e. they do not rapidly evolve away from their input values. In the 
future, however, experimental constraints for the nuclear gluons in this region of 
x could be obtained from the measurements of dileptons (originating from the 
decays of the D mesons) by the NA60 experiment at the SPS P7| . 



2.3 The EKS98 parametrization 

We have discussed above how to get constraints for the nPDF. In the DGLAP frame- 
work the problem boils down to determining the input distributions, i.e. the input 
ratios Rq(x,Qq), Ry(x,Q%) and R^(x,Qq). In practise, the extraction of the input 
ratios involves a recursive procedure: first the DGLAP evolution is performed with 



S 



some input distributions, then a comparison with the data is made at various values of 
x and Q 2 , after which the input distributions are changed in order to achieve a better 
agreement with the data. This iterative procedure is repeated until the best set of 
Rq(x, Qq) Ry(x,Ql) and R A (x, Ql) is found. The details can be found in 0. 

It is clear that the constraints always restrict the absolute nPDF ff' A (see e.g. Eqs. 
(|) and (||)). Constraints for the ratios Rf depend on the chosen set of parton distri- 
butions of the free proton, in terms of which the ratios in Eq. (|6]) are defined. Ideally 
of course there exists only one best set {/f }, but in practise, as some uncertainties still 
appear also in the free proton level, several sets are in use, (MRS, CTEQ, GRV, etc.). 
To test how strongly the obtained nuclear effects might depend on the choice of the 
parton distributions of the free proton, we have repeated the analysis of |7| in || by 
using the CTEQ4L |44j distributions instead of the GRVLO distributions [ 43j . In spite 



of the quite large differences in the gluon content (at small x there are fewer gluons 
in CTEQ4L) and in the sea quarks (more flavour asymmetry in CTEQ4L) between 
these sets, the differences between -R^ GRV and -R^ CTE q remained within a few percents 
only. Accepting this range of uncertainty, it is meaningful to prepare a "universal" 
parametrization for the ratios Rf(x,Q 2 ). This task was performed in ||, and the 
parametrization, called "EKS98" , is available for public use in the www [[45] and now 
also in the latest version of the CERN PDFLIB 



3 Comparison of different parametrizations 

Next we shall compare the EKS98-parametrization with two other parametrizations of 



the nuclear effects used in the literature, the default one in HIJING [|13j and the one 
prepared in one of the Hard Probe Workshops [Q. Especially, we shall focus on the 
comparison of all these against the data. 

3.1 The HIJING parametrization 

The HIJING parametrization of the nuclear modifications of the parton distributions 
given in Ref. [13| is the following: 



Rj 2 ( x ) = l+1.19(lnA) 1 / 6 [x 3 -1.5(x +x L )x 2 +3x x L x]-[a A - *\ A ■ ^ ^K^, 



\n{A + 1) 

(14) 

where a a = 0.1 (A 1 / 3 — 1), xq = 0.1 and xl = 0.7. The parametrization is based on fits 
to the old EMC data on Rp 2 |]18| . In particular, in the applications of this parametriza- 



tion, it is assumed that the modifications are identical for all parton flavours, and that 
the Q 2 dependence of the ratios is a negligible effect, Rf(x,Q 2 ) = Rp 2 (x). 
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3.2 The HPC parametrization 



An outcome of the Hard Probe Collaboration (HPC) meeting in ECT* at Trento (1995), 
was a parametrization [14] of Rp 2 (x). Similarly to the HIJING parametrization above, 
this parametrization assumes that the Q 2 dependence is a negligible effect. Possible 
differences between the modifications of different parton species are not considered, 
either. The HPC fit was motivated by a corresponding parametrization in and it 
was obtained by fitting the re-analysed SLAC data [pi 
The functional form of the HPC fit is 



and re-analyzed NMC data 



R a fM 



1 + c D c A (l/x - l/x sh ) 
Ksh l + c A A^(l/x-l/x sh ) 



~x, 



Pi 



X ^ 3^sh 



X s h < X < Xf 



Xf < X < 1 



(15) 



where the different regions are matched together by setting i? s h = a emc — & em c^sh and 
Rf — a cmc — KmcXf. In the EMC region a emc = 1 + fremc^emc- The A dependence of 
6„ m , is b cmc = p cmc [l - A- 1 / 3 - 1.145A" 2 / 3 + 0.93A" 1 + 0.88A" 4 / 3 - 0.59A" 5 / 3 ] from 
p8| . It was assumed (based on the data and the assumed Q 2 independence) that 
i emc ) = 1 and that the location of the EMC minimum is at Xf independently of A. 



J cmc 

Ref. 



The eight fit parameters obtain the following values: p A = 0.10011, c A = 0.0127343, 
0.154037, x cmc = 0.275097, p cmc = 0.525080, x f = 0.742059 and 



Pi z 



= 1.05570 x sh 
0.320992. 



3.3 The Comparison 

The results of the DGLAP analysis [§ for the ratios R&(x, Q 2 ), Rf(x, Q 2 ), R${x, Q 2 ), 
and R% at scales Q 2 = 2.25 GeV 2 (thin solid lines), 14.7 GeV 2 (dashed), 108 GeV 2 
(dotted-dashed) and 10000 GeV 2 (dashed) are compared with the HIJING parametriza- 
tion (thick solid) and with the HPC parametrization (thick dashed) in Fig. 2. This fig- 
ure is to illustrate two points: first, the Q 2 dependence becomes a non-negligible effect 
at small values of x. Secondly, the available experimental and sum-rule constraints lead 
to mutually quite different modifications for the valence quarks, sea quarks and gluons. 
For A = 208, the HIJING parametrization underestimates the baryon number sum rule 
by 18 % (24 %) and the momentum sum rule by 13 % (23 %) at Q 2 = Q 2 = 2.25 GeV 2 
(Q 2 = 10 4 GeV 2 ). For the HPC parametrization the corresponding figures are 5 % 
(12 %) for the baryon number and 8 % (12 %) for the momentum sum rule. The ef- 
fects of choosing different (lowest order) free proton parton distributions remain within 
one unit of percent. The deficit can easily be understood from the figure by comparing 
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Figure 2. The nuclear modifications of 
parton densities for a heavy nucleus A = 
208 (isoscalar). The ratios Rq(x,Q 2 ), 
R$(x,Q 2 ), R$(x,Q 2 ), and R%(x,Q 2 ) 
from Ref. @ are denoted by EKS98, 
(thin lines, plotted at the fixed values of 
Q 2 /GeV 2 indicated on the left). The Q 2 
independent parametrizations of the ratio 
Rj^ix) of HIJING (thick solid line) and 
that of HPC (thick dashed line) remain 
unchanged from panel to panel. The no- 
tation is the same in all panels. 



with the EKS98 in which these conservation laws are met within an accuracy of ~ 1 % 
at all Q 2 for any (lowest order) PDF set of the free proton. 

Fig. 3 shows the comparison of the EKS98 (open squares), HIJING (open diamonds) 
and HPC (open triangles) parametrizations with the NMC data for Rp 2 (filled 
circles). The EKS98 results are computed by using CTEQ5L set of parton distributions 
|30| , |46[| and at scales Q 2 corresponding to the (Q 2 ) measured for each x (see Fig. 1). 
The differences remain small between the EKS98, HPC and the data. HIJING agrees 
with the data for carbon but overestimates the A-dependence of shadowing for other 
nuclei. The error bars in the data represent the statistical and systematic errors added 
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Figure 3. Comparison of the EKS98 
(open squares), HIJING (open diamonds) 
and HPC (open triangles) parametriza- 
tions with the NMC data (filled circles 
with error bars) for F^/F® for helium, 
carbon and calcium (2*|] . For EKS98 the 
comparison is only made for the experi- 
mental points above Q 2 = 2.25 GeV 2 . No- 
tice that the vertical scales are different in 
each panel. 



5 10" 



5 10- 



5 10" 



X 



in quadrature. 

Next, in Fig. 4 we show the comparison of EKS98, HIJING and HPC with the 
A systematics of Rp 2 measured by the NMC [p5|l . Again, the data is shown by the 
filled circles with statistical errors (the inner error bars) and with the statistical and 
systematic errors added in quadrature (the outer error bars). The notation is the same 
as in the previous figure, and the EKS98 results are again computed at the scales 
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Figure 4. Comparison of the EKS98 (open squares), HIJING (open diamonds) and HPC (open 
triangles) parametrizations with the NMC data (filled circles with error bars) for F^/F^ 



Q 2 = (Q 2 ) for each value of x given by the experiment. This figure shows that the 
A dependence of shadowing is clearly too strong in HIJING. It should also be noted 
that the data of Fig. 4 was not yet available for the HPC fit but that the HPC 
parametrization falls nevertheless fairly close to the data. Notice that of the data sets 
shown in Fig. 4, the one for Sn/C gives the most stringent constraint for any fit, due 
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to the smallest error bars. 

In Fig. 5, we have plotted the Q 2 dependence of the ratio F^/ -F 2 C at fixed values of x 
corresponding to those in the NMC data |f2B] . The EKS98 (with CTEQ5L distributions) 
is shown by the solid lines, HIJING by the dotted lines and HPC with the dashed lines. 
The data is shown by the open squares with (statistical) error bars. At small values of 
x the Q 2 dependence is not a negligible effect. 
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Figure 5. The Q dependence of the ratio F2/F2 for fixed values of x. The EKS98 results 
are the solid lines, HIJING parametrization the dotted line and the HPC parametrizaton the 
dashed line. The NMC data [26] is shown by the filled circles with (statistical) error bars. 



Finally, in Fig. 6 we show the comparison of EKS98, HIJING and HPC with the 
E772 DY data 



33 1 in pA collisions. In the computation of the ratio from Eq. 
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Figure 6. The Drell-Yan ratio R J ^ )Y { X 2-, (Q 2 )) computed from Eq. (||). The E772 data [ 33 1 are 
shown by the filled circles with error bars, the EKS98 results by the open squares, the HIJING 
parametrization by the open diamonds, and the HPC parametrization by open triangles. 



the CTEQ5L distributions [0, [|6| of the free proton have been used. The scales Q 2 
for each x 2 used in computing the EKS98 results are those in Fig. 1. Also here the 
conclusion is that the HIJING parametrization clearly overestimates the A dependence 
of shadowing. 



4 Conclusions 

In this note we have discussed the determination of the nuclear parton distributions in 
the lowest order leading twist DGLAP framework. We have shown which kinematical 
range and which combinations of the parton distributions are probed in the measure- 
ments of the deeply inelastic I A scatterings and the Drell-Yan process in pA collisions, 
and to what extent the measurements offer constraints for the input distributions. The 
remaining uncertainties are also discussed. 

To demonstrate the differences between the different parametrizations of the nu- 
clear effects used in the literature, we have compared the EKS98, HIJING and HPC 
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parametrizations with each other and, most importantly, with the data. The HIJING 
parametrization for Rp 2 clearly overestimates the A dependence of nuclear shadowing, 
and leads to a contradiction with the data. Our conclusion therefore is that for de- 
tailed studies of shadowing effects the HIJING parametrization ]13j is not an adequate 
description of the nuclear modifications of quark and antiquark distributions. If the 
flavour and scale dependence of the nuclear effects could be neglected for some applica- 
tions, we note that the HPC parametrization gives a better representation of the data 
for the ratios Rp 2 than the one in For such a use the HPC fit should, however, be 
redone by including the extensive A systematics provided by the NMC measurements 
25|j : especially the data on the ratio F 2 Sn /F 2 c offers an additional constraint for the fit. 



In the leading twist framework also the conservation of baryon number and mo- 
mentum should be required, since all the nuclear effects are contained in the parton 
distributions. If so done, the modifications for the valence quarks, sea quarks and glu- 
ons differ from each other, as is demonstrated in the DGLAP analyses P, H, 0, |, Hi- 
If the HIJING or HPC parametrizations of Rp 2 (x) are directly used to modify the dis- 
tribution of all parton flavors, the baryon number will be underestimated for A = 208 
by 18.. .24 % with HIJING and 5.. .12 % with HPC at scales Q 2 = 2.25.. .10 4 GeV 2 . 

The NMC measurements J26| have revealed a clear Q 2 dependence in the ratio 
F 2 Sn /F 2 c , at x ~ 0.01. Naturally, scale independent parametrizations cannot repro- 
duce the observed behaviour. In a DGLAP framework, the Q 2 dependence of the ratio 
F 2 Sn /F 2 c can be used to constrain the nuclear gluon distributions, as first suggested 
in Ref. [ji^] . In the DGLAP analysis it was shown that the constraints obtained 
for the gluons [J2"| are also consistent with the momentum sum rule, i.e. that a fairly 
strong antishadowing of gluons appears at the input scale Ql ~ 2 GeV 2 . The HIJING 
parametrization underestimates the momentum sum for A = 208 by 13. ..23 % and the 
HPC by 8. ..12 % at scales Q 2 = 2.25.. .10 4 GeV 2 . In addition, within the DGLAP 



analysis, the data |26| seems to rule out the case such as the Ansatz 2 in Ref. Q, 
where gluons would have clearly stronger shadowing than that observed in Rp. 

Finally, the analysis of Ref. @ which lead to the EKS98 parametrization || should be 
improved within the DGLAP framework in obvious ways: the fitting procedure should 
be automated for inclusion of new data sets in the future, the parameter space should 
be more thoroughly explored to estimate better the uncertainties (see also ||) , and 
next-to- leading order analysis should be performed. We do not, however, expect the 
results change very much from the EKS98 j7|, |§. One should also consider expansions 



of the DGLAP framework, such as an inclusion of the recombination terms [|{| [5(J in 
the evolution equations @, ||, parton saturation phenomena at small values of x [49]- 
53 ], and possible higher-twist effects in the cross sections [ |TT| [TTJ, especially at lower 



scales. Also input from the different models for the origin of the nuclear effects can be 
considered. However, the comparison with the data should remain as the key feature 
of the analysis, since it is (at least presently) not possible to compute the absolute 
nuclear parton distributions from first principles. 
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